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CN . Abstract 

m 

• I First measurements are presented of the diffractive cross section aep^eXY at centre-of- 

Q ■ mass energies ^/s of 225 and 252 GeV, together with a precise new measurement at ^/s 

of 319 GeV, using data taken with the HI detector in the years 2006 and 2007. Together 
with previous HI data at ^/s of 301 GeV, the measurements are used to extract the diffrac- 
tive longitudinal structure function F^^ in the range of photon virtualities 4.0 < Q^ < 
^ . 44.0 GeV^ and fractional proton longitudinal momentum loss 5 • 10~^ < xp < 3 • 10^^. 

H I The measured F[' is compared with leading twist predictions based on diffractive parton 

densities extracted in NLO QCD fits to previous measurements of diffractive Deep-Inelastic 
Scattering and with a model which additionally includes a higher twist contribution derived 
from a colour dipole approach. The ratio of the diffractive cross section induced by lon- 
gitudinally polarised photons to that for transversely polarised photons is extracted and 
compared with the analogous quantity for inclusive Deep-Inelastic Scattering 
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1 Introduction 

The observation that a significant subset of Deep-Inelastic Scattering (DIS) events at HERA 
contain a large gap in activity in the forward region [IJ prompted much theoretical and experi- 
mental work. Such large rapidity gap topologies signify a colour singlet or diffractive exchange 
and HERA has proved to be a rich environment for their study. In particular, the study of 
diffractive DIS (DDIS), both inclusive and exclusive, has supplied a wealth of experimental 
data with a hard scale given by the photon virtuality, stimulating the theoretical understanding 
of diffraction in terms of perturbative quantum chromodynamics (QCD). 

It has been shown that the neutral current DDIS process ep — )■ eXp at HERA obeys a QCD 
factorisation theorem [2J. This allows for a description of DDIS in terms of parton densities con- 
voluted with hard scattering matrix elements. The diffractive parton density functions (DPDFs) 
depend on four kinematic variables, so an additional assumption is often made whereby the pro- 
ton vertex dynamics factorise from the vertex of the hard scattering, as shown in figured] While 
this proton vertex factorisation has no complete foundation in theory, measurements of DDIS 
from both HI [|3]-|5l and ZEUS [6] show that it holds well enough such that next-to-leading or- 
der (NLO) QCD fits can be made to the data [[3ll7]-l9l|. The DPDFs then depend only on the scale 
Q"^ and the fraction z of the total longitudinal momentum of the diffractive exchange which is 
carried by the parton entering the hard scattering. 

Measurements of the dijet cross section in DDIS allow tests of the DPDFs extracted in fits 
to inclusive DDIS data. This process, which is known to be dominated by boson-gluon fusion, 
is particularly sensitive to the poorly known gluon DPDF at large z and has thus been used 
successfully to distinguish between different DPDF sets [9|. DDIS events containing charm 
particles in the final state have similarly been used to test the gluon DPDF [fTOl. 

As in the inclusive DIS case, the cross section for DDIS can be expressed in terms of a linear 
combination of structure functions, F^ and F^ [fTTI . While F^ describes the total photon- 
proton process, F^ is only sensitive to the longitudinally polarised photon contribution. As for 
its inclusive counterpart, F^ is thus zero in the quark-parton model, but may acquire a non-zero 
value, < F[' < F^ in QCD, with leading twist contributions dependent on both the diffractive 
quark and gluon densities [12J. A measurement of Ff provides a powerful independent tool to 
verify our understanding of the underlying dynamics of diffraction up to NLO in QCD and to 
test the DPDFs. This is particularly important at the lowest z values, where direct information 
on the gluon density cannot be obtained from dijet data due to kinematic limitations and where 
novel effects such as parton saturation [frBil or non-DGLAP dynamics I114[|15ll are most likely to 
become important. 

Previous attempts to measure F^ [IS [HI have exploited the azimuthal decorrelation be- 
tween the proton and electron scattering planes expected due to interference between the am- 
plitudes for transverse and longitudinal photon polarisations [[T71 . However, due to the rel- 
atively poor statistical precision of the measurement, the results were consistent with zero. 
The HI collaboration has recently published measurements of the inclusive structure function 
Fl{x, Q"^) [I18U19 1 using the centre-of-mass energy dependence of the DIS cross section at fixed 
X and Q^. A similar approach has been proposed to extract F^ [|20l . 

In addition to measuring Ff itself, it is interesting to compare the relative sizes of the 
diffractive cross sections induced by transversely and longitudinally polarised virtual photons. 
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Figure 1: A diagram of the diffractive DIS process ep — >■ eXp or ep — )■ eXY. The dotted line 
indicates where the diagram can be divided under the assumption of proton vertex factorisation. 



This comparison has previously been made for inclusive DIS and exclusive vector meson pro- 
duction through the study of the photoabsorption ratio, R = o-l/ctt , where ai and ar are 
the cross sections for the scattering of longitudinally and transversely polarised photons, re- 
spectively. Whilst R is only weakly dependent on kinematic variables in the DIS regime for 
inclusive cross sections lfT8ll2Tll . a strong dependence on Q^ is observed for vector meson pro- 
duction [22], the longitudinally polarised photon cross section becoming much larger than its 
transverse counterpart at large Q^. Since DDIS incorporates vector meson production and re- 
lated processes at large z, but exhibits kinematic dependences which are similar to those of 
inclusive DIS at low z, it is not easy to predict its photoabsorption ratio. By analogy with the 
inclusive DIS case, we define R^ = F^ /(F,f — F£) for diffraction. The double ratio R^ /R 
thus measures the relative importance of the longitudinally and transversely polarised photon 
cross sections in diffractive compared with inclusive scattering. 

In this analysis, positron-proton collision data taken at different proton beam energies with 
the HI detector at HERA in the years 2006 and 2007 are used to measure the diffractive cross 
section at intermediate and large inelasticities y. Dedicated low and medium energy (LME) 
data with proton beam energies of Ep = 460 and 575 GeV are analysed together with data 
at the nominal beam energy of 920 GeV. Previously published data at a proton beam energy 
Ep = 820 GeV [3] are used in addition. The positron beam energy is 27.6 GeV in all cases. 
These cross sections are used to extract F^ together with the ratio R^ and the double ratio 
R^'/R. 



2 Kinematics and cross section definition 



The kinematic variables used to describe inclusive DIS are the virtuality of the exchanged boson 
Q^, the Bjorken scaling variable x and the inelasticity variable y, defined as: 

Q'' = -q^ = -(k-k'f x = -^ y = ^-^, (1) 

^ ^ ^ ' 2P-q ^ P-k' ^ 

where k and k' are the four-momenta of the incoming and outgoing positrons, respectively, 
and P is the four-momentum of the incoming proton. They are related to s, the square of the 
centre-of-mass energy, by Q"^ = sxy. 

In diffractive events, the hadronic final state can be divided into two systems X and Y 
which are separated by the largest gap in rapidity. A diagram for the DDIS process is shown 
in figure [TJ The system Y is either the elastically scattered proton, which is the dominant state 
in the kinematic range studied here, or its low mass excitations. In addition to the standard 
DIS variables and the squared four-momentum transfer at the proton vertex, t, the kinematic 
variables xp and /3 are useful in describing the diffractive DIS interaction. They are defined as: 

q-jP-py) r, Q" ... 

q- P 2q- {P-Py) 

where py is the four momentum of the elastically scattered proton or of its low mass excitation. 
The variable xjp is the longitudinal momentum fraction of the proton carried by the diffractive 
exchange and (3 is the longitudinal momentum fraction of the struck quark with respect to the 
diffractive exchange, such that x = X]p(3. In the simple quark-parton model, /3 = z, while for 
higher order processes, < /3 < z. The results are discussed in terms of a diffractive reduced 
cross section, (t^{(3, Q^, xp), related to the measured differential cross section by: 
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^ ■ F+ ■ a^^^\xjp,P,Q^), (3) 



dxjpd/3dQ2 /3Q4 

where F+ = 1 + (1 — y)^. The diffractive reduced cross section is related to the diffractive 
structure functions by: 

a,^(3)(x^,/3,Q2) = F^^'\xjp,P,Q') - ^F^^'\xjp,P,Q^). (4) 

Due to the suppression term y'^/Yj^, the diffractive reduced cross section is only sensitive to F^ 
at large values of y. 

As the final state system Y is not measured in this analysis, the cross section is integrated 
over ranges in its mass My and in t. These ranges are chosen to be 

My < 1.6 GeV, |t| < 1.0 GeV^ (5) 

corresponding to the acceptance of the HI detector in the forward direction and for consistency 
with previous measurements. 



3 Models of Ff 



The relationships between the diffractive structure functions and the DPDFs have been shown 
to be analogous to those of the inclusive case in the limit where the proton mass and t may 
be neglected compared with other relevant scales in the interaction [TT] . The diffractive DIS 
structure function F^ is then directly sensitive to the singlet quark DPDF and the scaling vi- 
olations, dF^ /d In Q"^, provide a measure of the gluon DPDF. NLO QCD fits to a^ at low to 
intermediate y values, sometimes supplemented by dijet data, thus provide DPDFs which lead 
to predictions of F^ at leading twist. By analogy with the inclusive case [[T2l|23l and assuming 
coUinear factorisation ^, the NLO expression for F^ in the MS scheme is 



?D{3) 



{fi,Q',Xjp) 



27r 



dz 



^ E^^^(f'^'^^) + /^(f'^'^^)(i- 



k={q,q} 



, (6) 



where fg and fg are the quark and gluon DPDFs and e^ is the electric charge of quark flavour 
k. At the relatively large (3 values at which Ff can be measured at HERA, both the quark and 
the gluon densities are predicted to make important contributions to F^ , despite the dominant 
role played by gluons in DDIS in general |i3l|7|. 

In this paper, the Ff measurement is compared with predictions derived from two NLO 
QCD fits to inclusive DDIS a^ data [3J, which are labelled 'HI 2006 DPDF Fit A and 'HI 
2006 DPDF Fit B'. Proton vertex factorisation is assumed in both cases and the diffractive 
quark densities are very similar in the two fits. However, the two DPDF fits differ in their pa- 
rameterisations of the gluon density, which leads to considerable differences at large fractional 
momenta z [|3||, where the constraints from inclusive DDIS data are poor. Corresponding differ- 
ences are visible between the Fit A and Fit B predictions for Ff . The 'HI 2006 DPDF Fit B' 
DPDFs give the better description of diffractive dijet production at HERA [9] and are therefore 
used as the default here. 



A complementary approach to modelling diffractive DIS is offered by dipole models 11241 
|25^. Viewed in the proton rest frame, the incoming virtual photon fluctuates into a qq pair or 
higher multiplicity state, whose scattering strength from the target is governed by a universal 
dipole cross section. Dipole models which are applicable to DDIS generally contain three 
contributions II24[|26II : leading twist terms corresponding to the scattering of qq and qqg dipoles 
derived from fluctuations of transversely polarised photons, and a higher twist contribution 
(suppressed like l/Q^) in which qq dipoles are obtained from longitudinally polarised photons. 
Dipole models thus tend to neglect the leading twist contribution to Ff which emerges naturally 
from NLO DPDF fits. However, the higher twist contribution to Ff is of particular interest, 
since it can be predicted in perturbative QCD [|27|| , by coupling a qq dipole to a two-gluon 
exchange in a similar phenomenology to that successfully applied to vector meson cross sections 
at HERA [28] . In many dipole-inspired models, this higher twist component is the dominant 
feature of a^ at large (3 and low-to-moderate Q"^. 

In a recent hybrid approach to fitting a^ [29] (labelled 'Golec-Biernat & Luszczak' here), 
the leading and higher twist contributions to F^ are included simultaneously. A parametrisation 
similar to that in jSl is used for the diffractive quark and gluon DPDFs, but the higher twist 



longitudinal photon contribution is also included via the parametrisation employed in [24] . The 
quality of the fit to the a^ data is similar with and without the higher twist term. However, its 
inclusion leads to a sizeable effect on the diffractive gluon density at large fractional momenta 
and the higher twist contribution dominates the resulting predictions for Ff for (3^0.6 at the 
lowest Q'^ values considered here. 



4 Experimental Method 
4.1 HI detector 

A detailed description of the HI detector can be found elsewhere [30] and only the components 
essential to the present analysis are briefly described here. The origin of the HI coordinate 
system is the nominal ep interaction point at the centre of the detector, with the direction of 
the proton beam defining the positive 2;-axis (forward direction). The polar angle (9) is defined 
with respect to this axis and the pseudorapidity is defined as 77 = — In tan{9/2). The azimuthal 
angle defines the particle direction in the transverse plane. 

The analysis uses several of the tracking detectors of HI, relying primarily on the two con- 
centric central jet chambers (CJC) and the central silicon tracker (CST) [|3TI . which measure the 
transverse momenta of charged particles in the angular range 20° < 9 < 160°, together with the 
backward silicon tracker (BST), which is positioned around the beam-pipe in the backward di- 
rection. Complementary tracking information is obtained from the z drift chamber COZ, which 
is located in between the two cylinders of the CJC, the forward silicon tracker (FST) and the 
forward tracking detector (FTD). The central inner proportional chamber (CIP) [32] provides 
trigger information on central tracks, the FST and BST are used to improve the overall vertex 
reconstruction and the FTD is used to improve the hadronic final state reconstruction of low 
momentum particles in the forward direction. 

In the backward region —4.0 < r] < —1.4, a lead- scintillating fibre calorimeter (SpaCal) is 
used for the identification and measurement of the scattered positron, with an energy resolution 
for electromagnetic showers of cr(ii^)/ii^ ~ 7 .1% / y^ E / GeV (Bl%. Importantly, it also provides 
a trigger down to positron energies of 2 GeV. The hadronic section of the SpaCal is used 
in the reconstruction of the hadronic final state, especially at the high y values accessed in 
this analysis. The liquid argon (LAr) calorimeter covers the range —1.5 < 77 < 3.4 and is 
also used in this analysis in the reconstruction of the hadronic final state. It has an energy 
resolution of a{E)/E ~ 50%/y/EjGeV for hadronic showers, as obtained from test beam 
measurements [|33l . 

Several of the forward detectors of HI are used in conjunction with the LAr to determine 
whether or not an event contains a large rapidity gap close to the outgoing proton direction. 
The forward muon detector (FMD) comprises two sets of three drift chambers, separated by a 
toroidal magnet, covering the range 1.9 < r/ < 3.7. Only the three layers closest to the interac- 
tion region are considered in this analysis. A dedicated reconstruction algorithm efficiently de- 
tects secondary particles produced through the interactions of proton dissociation products with 
the beam-pipe or other accelerator elements, giving the FMD an effective coverage extending 
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to around rj = 6.5. The Plug is a calorimeter consisting of four double layers of scintillator and 
lead absorber, read out by photomultipliers. It is situated at z = 4.9 m and covers the range 
3.5 < T] < 5.5. The final forward detector component used in the analysis is one station of 
the forward tagging system (FTS), consisting of scintillators situated around the beam-pipe at 
z = 28ia covering approximately 6.0 < i] < 7.5. 

Positrons scattered through very small polar angles can be detected with a calorimeter 
(ETAG) placed at 2; = — 6 m downstream in the positron beam direction. The luminosity 
is determined from the Bethe-Heitler scattering process, which is measured using a photon 
calorimeter at z = —103 m. 



4.2 Data samples 

Three samples are analysed to provide data at different centre-of-mass energies in different 
kinematic ranges, as shown in table [1] 



(GeV) 


(GeV) 


Q^ range 

(GeV^) 


y range 


Luminosity 

(pb'^) 


460 
575 
920 


225 
252 
319 


2.5 < Q" < 100 
2.5 < Q2 < 100 
7.0 < Q2 < 100 


0.1 <y<0.9 
0.1 <y<0.9 
0.1 <y<0.56 


8.5 

5.2 

126.8 



Table 1 : Summary of the data samples used in the analysis. 

In addition to these data, cross section measurements at Ep = 820 GeV from a previous HI 
publication [[3l are used to extract Ff^ in the same kinematic range. 



4.3 Event selection 



Dedicated 'high y' triggers are used for the LME datasets in order to allow triggering on energy 
depositions as low as 2 GeV in the SpaCal. For y > 0.6 (0.56) in the 460 (575) GeV data, the 
SpaCal trigger decision is combined with information from the BST or CIP in order to reduce 
the rate. For lower y values, corresponding to high energy depositions in the SpaCal, triggers 
based on SpaCal-only information are used for all three datasets. The combined efficiency of 
the LME high y triggers is around 99% for positron energies above 3 GeV, as monitored with 
independent triggers. The data are corrected for this inefficiency, which has a small depen- 
dence on the radial position of the scattered positron in the SpaCal, Rspacah due to the track 
requirement. The combination of SpaCal-only triggers used has a negligibly small inefficiency. 

The event selection is based on the identification of the scattered positron as a localised 
energy deposition, a cluster, of more than 3.4(12.0) GeV in the SpaCal in the LME (920 GeV) 
data. Backgrounds due predominantly to photoproduction processes, where the scattered positron 
is lost down the beam-pipe, are reduced by requiring that the logarithmic energy-weighted clus- 
ter radius, riog, is smaller than 5 cm and that the energy measured in the hadronic section of the 
SpaCal associated with the cluster is less than 15% of the cluster energy. If the highest energy 



cluster fails to fulfill these selection criteria, the second and third highest energy clusters are 
considered in turn. QED Compton contributions, ep — )■ e^yp, are suppressed by rejecting events 
with two back-to-back clusters. 

For the LME data, the background is further reduced by demanding a 'linked track' that 
can be extrapolated to the SpaCal cluster within a radial distance of 3 cm. The linked track is 
reconstructed using a dedicated algorithm incorporating information from both the CJC and the 
BST [34]. Geometrical cuts are applied to keep the tracking acceptance high and track quality 
requirements are applied, reflecting the geometry of these detectors. 

In order to further reject background, a reconstructed event vertex is required to lie within 
35 cm of the nominal interaction point for all data samples. In order to guarantee a high vertex- 
finding efficiency, the measurement is restricted to the kinematic range y > 0.1. An algo- 
rithm combining calorimeter and tracking information, which optimises precision while avoid- 
ing double-counting, is used to reconstruct the four vector of the hadronic final state (HFS) 
particles [35|. For all datasets, the quantity Sj(i? — pz)i, where the sum is over the energy E 
minus the longitudinal momentum p^ of all final state particles including the scattered positron, 
is required to be greater than 35 GeV. This quantity should peak at twice the incident positron 
energy, i.e. 55 GeV, for fully reconstructed DIS and DDIS events alike. This completes the 
background rejection criteria of the inclusive event selection. 

At low positron energies, the photoproduction background remains large after all cuts. Fol- 
lowing the procedure explained in [|18il . this residual background is estimated from the number 
of events N^c passing the full analysis selection and having a negatively charged track linked 
to the SpaCal cluster. The photoproduction background is expected to be approximately charge 
symmetric and therefore corresponds to approximately 2N]yc- However, a small asymmetry in 
its charge composition has previously been measured ifTSl . Thus the photoproduction estimate 
is 1.98Nwc, which is statistically subtracted from the sample. 

Diffractive DIS events are selected as a subsample of the inclusive DIS event sample on 
the basis of a large rapidity gap in the forward direction. The pseudorapidity r]max of the 
forward-most energy deposit above 800 MeV in the LAr calorimeter is required to be less 
than 3.3. In addition, the FMD, Plug and FTS are required to have no discernible signal above 
their typical noise levels. The combined efficiency for rejecting proton dissociative events with 
My ^1.6 GeV is greater than 99%. These requirements select a subsample of events where the 
hadronic final state is separated into two systems X and F by a large rapidity gap. The system 
Y, which is predominantly a single proton, escapes undetected down the beam-pipe, while the 
system X is fully contained in the main HI detector. 

In order to maintain a high efficiency for the vertex reconstruction of the DDIS event sam- 
ple, an additional fiducial cut is required to avoid cases where both the final state system X and 
the positron are outside the acceptance of the CJC. The region where both Rgpacai < 40 cm 
and r]max < —1-7 is removed from the analysis, after which the vertex-efficiency is high and 
well understood throughout the measured phase space. Finally, there must be at least one re- 
constructed HFS particle to define the system X. 

The inclusive DIS event kinematics are reconstructed using different methods depending 
on the y range of a given dataset. For the LME data, only information from the reconstructed 
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scattered positron is used, as this method has the best resolution at large y: 

y = 1 - TT sm — Q = x = — . (7) 

Ee 2 1-y sy 

Here, E^ is the energy of the incident positron and E^ and 6^ are the energy and polar angle of 
the scattered positron, respectively. For the 920 GeV data, a method with better performance at 
low y is used [|36]| : 



.2,., n ., ^ r.2_ ^E!{l-y) ^ ^ Q^ 

tan^ 6e/2 sy 



y = yl + y,{l-y,) Q^ = ^^^1 3; = ^, (8) 



where yd = tan (7/2)/[tan {9e/2) + tan (7/2)] and 7 is the polar angle of the hadronic final 
state. 

The four momentum of the final state system X is reconstructed as the vector sum of all 
HFS particles. Its mass Mx is reconstructed as: 



Mx = f{Vma.)J{E' -Pi -Pi- PDu.S-. (9) 

V Vh 

where {E,px,Py,Pz)„ps denotes the four vector of the HFS and yh = {E — pz)jjpg/2Ee. The 
term y/yn improves the resolution and the function f{r]max) is determined from simulation and 
corrects for detector losses. The diffractive variables are then reconstructed as: 



4.4 Corrections to the data and simulations 

Monte Carlo (MC) simulations are used to correct the data for the detector effects of acceptance, 
inefficiencies, and migrations between measurement intervals. The DDIS signal is modelled 
for xjp < 0.15 using the RAPGAP [[371 generator, with HI 2006 DPDF Fit B [H as the input 
DPDFs. Higher order QCD radiation is modelled using initial and final state parton showers 
in the leading log{Q'^) approximation ['381. Hadronisation is simulated using the Lund string 
model |l39l as implemented in PYTHIA [40]. As RAPGAP is a leading order MC generator 
simulating only F^, the effect of Ff has been simulated by weighting RAPGAP events by the 
ratio aj^ /F^ as given at NLO by HI 2006 DPDF Fit B. This is important at high y in order to 
describe the data. At low Q'^, HI 2006 DPDF Fit B undershoots the data, as observed previously 
p] . RAPGAP is therefore reweighted for Q^ < 7 GeV^ by a parametrisation of the ratio of the 
previous data to HI 2006 DPDF Fit B. Resonant contributions to the diffractive cross section, 
important at low Q"^ and low Mx < 5 GeV, are modelled using the DIFFVM [HTI generator. 
The DIFFVM generator is also used to simulate proton dissociative events with My < 5 GeV 
to correct the measurements to the My and t ranges given in equation [5] under the assumption 
of proton vertex factorisation. The small non-diffractive DIS background from xp > 0.15 or 
My > 5 GeV is modelled using DJANGO JH, while the COMPTON program [431 is used to 
model the QED Compton process, important at very low Mx- 
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The generated events are passed through a full GEANT [|44l simulation of the HI detector. 
The simulated events are subjected to the same reconstruction and analysis chain as the data. 
More details of the analysis can be found in [i45l . 

Figure [2] shows the energy distributions for positron candidates in the LME datasets. In 
addition to the simulation described above, the photoproduction estimate using the number of 
candidates with the wrong charge, and the total background expectation are also shown. The 
data are well described down to positron energies of 3.4 GeV. 

The quality of the calibration of the system X, in the sensitive region at high y, is illustrated 
in figure|3l where Sj {E—pz)i peaks at the expected value of 55 GeV and is well described by the 
simulation. At large y, the hadronic energy measurement is strongly influenced by the hadronic 
energy response of the SpaCal, which has been calibrated using inclusive DIS events [|45l . The 
influence of varying the SpaCal hadronic energy scale by ±5% is indicated in the figure. 

The y, /3 and log(xjp) distributions in the data are compared with the total expectation in 
figure |4] for all three datasets. Again, the photoproduction estimate and the sum of all other 
background sources are also shown. The quality of the description is good in all cases. 



4.5 Cross section extraction 

The data are analysed in two Q^ ranges. For Q^ > 7 GeV^, data are available from all three 
datasets at Ep = 460, 575 and 920 GeV. For 2.5 < Q^ < 7.0 GeV^ only data from the 460 
and 575 GeV datasets are analysed. Previous measurements at Ep = 820 GeV jSl are used in 
addition in the Q^ and xp range of the LME data. The Q^, xp and (3 values of these published 
data have been adjusted to the values of the current analysis using a parameterisation of a^ 
derived from HI 2006 Fit B, a procedure which results in a systematic uncertainty of 1% at 
X]p = 0.003 and 3% at xp = 0.0005. The reduced cross section is extracted as a function of 
/3, Q^ and xp from measurements of the differential cross section according to equation |3l The 
Q^ and xp measurement intervals are large and have been optimised for the extraction of Ff 
in as broad a kinematic range as possible. 

The data are corrected for efficiencies and migrations between measurement intervals using 
the MC simulation described in section |4~4l The acceptance, as calculated from the MC model, 
is required to be above 20% for all points and is much larger than this except at the lowest Q^ 
and Xp. Purity and stabilitjo are larger than 50% in all bins. For the LME data, the estimate 
of the photoproduction background using the number of candidates with the wrong charge, 
Nwc, is subtracted bin-by-bin for y > 0.6, while below this value the background is negligible. 
Inclusive DIS and QED-Compton contributions are also subtracted bin-by-bin using the MC 
simulations described in section |44l The parametrisation of a^ using HI 2006 DPDF Fit B is 
used to correct the data to the central Q^, xp and /3 values quoted. As /3 — )■ 1, the shape of the 
cross section is largely unconstrained by data and varies quickly due to resonant contributions, 
making the correction to a single point in the phase space problematic. Thus, for /3 > 0.9, the 
average cross section in that interval is given. 



' Purity is defined as the fraction of reconstructed MC events in a measurement interval which also originated 
in the same interval at the hadron level. Stability is the fraction of MC events in a measurement interval at the 
hadron level which are also reconstructed in that interval. 
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The diffractive reduced cross section is integrated over the My and t ranges given in equa- 
tion[5l DIFFVM is used to calculate the correction to this phase space, which varies with proton 
beam energy. The correction factors are 1.04, 1.06 and 1.15 for the 460, 575 and 920 GeV data, 
respectively. 

For use in forming the ratio R^ /R, inclusive cross sections are measured in the same bin- 
ning scheme as is used for the diffractive measurement, using the procedure described in IJTSl . 
As the statistics for the inclusive DIS sample are larger, the background subtraction is more so- 
phisticated. The number A^^ of events passing the full analysis selection and having a signal in 
the ETAG photoproduction tagger and a negatively charged linked track associated to a SpaCal 
cluster provides another estimate of the photoproduction background. For the 460(575) GeV 
data, at low y < 0.6(0.56), the photoproduction estimate uses Nt, while for higher y the pho- 
toproduction background is estimated using the number of candidates with the wrong charge, 
Nwc- For the 920 GeV data, the estimate based on positron-tagged events is used for all y. 



4.6 Systematic uncertainties 

A full systematic error analysis is performed, which carefully considers correlations between 
measurement intervals and data at different centre-of-mass energies. The sources of systematic 
uncertainty that have correlations between cross section measurement points at different Ep 
values are as follows. 

• The uncertainty on the electromagnetic energy scale of the SpaCal is 0.2% at the kine- 
matic peak of _Eg = 27.6 GeV, increasing linearly such that it would be 1% at E'^ = 

IGeV. 

• The possible bias in 9e is estimated using the mean difference in polar angle between the 
linked track and the SpaCal cluster, which is measured to be less than 1 mrad. 

• Noise is simulated in the LAr calorimeter using randomly-triggered events. The fraction 
of energy identified and subtracted as noise is known to a precision of 15%. 

• The hadronic section of the SpaCal is calibrated to a precision of 5%. The uncertainty on 
the hadronic energy scale of the LAr calorimeter is 2% and is found to have only a small 
effect on the cross sections in the present analysis. 

• The efficiency of the cut on the logarithmic energy-weighted cluster radius, r^og, is known 
toaprecisionofO.5%, 1.5%and3%for0.6 <y< 0.7,0.7 <y < 0.8 and 0.8 <y< 0.9, 
respectively. 

• The charge asymmetry in the lepton candidates from photoproduction background events 
of 0.98 is known to 4% precision filSll . 



D 
L 



The RAPGAP MC is weighted by the ratio of a^ /F.p in order to describe the data at 
high y. The associated uncertainty is evaluated by replacing Ff in the expression used 
for a^ in the reweighting procedure (equation S]) by either 0.5 ■ F^ or 1.5 ■ Fj 
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The kinematic dependences of the model used to correct the data are generally well con- 
strained from previous measurements. The uncertainties on the t, (3 and xjp dependences 
are evaluated by weighting the generator-level kinematics by e^\ /3='=°°^, (1 — /3)^°°^ and 
(l/xjp)^^^. The effects of weighting in t and (1 — /?) are found to have a negligible effect 
on the measured cross sections. 

The uncertainty due to the resonant contributions modelled by DIFFVM is evaluated by 
calculating the change in acceptance when including this contribution in the simulation 
or not. 

The non-diffractive DIS and QED-Compton backgrounds are modelled using MC simu- 
lations and are statistically subtracted from the data. The non-diffractive DIS background 
has a negligible effect in this analysis except at the highest xp. The QED-Compton events 
are only relevant for Mx -> 0. The normalisations of these backgrounds are controlled 
at the level of 100% and 30%, respectively. 

The corrections due to the finite measurement intervals (bin-centre corrections) are sub- 
ject to an uncertainty, which is evaluated from the change in these corrections when this 
procedure is carried out using the HI 2006 DPDF Fit A and Fit B parameterisations of the 
reduced cross section. The uncertainty is very small except at large (3, where the shape of 
a^ is not well constrained, and at low f3, corresponding to high y. 



Sources of experimental uncertainty which lead to systematic errors which are not correlated 
ween dc 
following. 



between data at different Ep values are the statistical errors of the MC simulations and the 



• The vertex reconstruction efficiency of the CJC is controlled to the level of 2% for xjp > 

10-=^ and 10% for lO^^ < xp < lO^l 

• The trigger efficiency is Z 99% and measured with a precision of 1% using independently- 
triggered data. 

• The uncertainty in the efficiency of linking a track to a SpaCal cluster is 1.5%. 

• The uncertainty on the efficiency of the forward detector selection for rejecting proton 
dissociative events is 0.5% [45 J. 

The model dependent uncertainties on the factors applied in correcting the measurements to 
the My and t ranges given in equation[5]are evaluated using the method described in [[3l. The re- 
sulting normalisation uncertainities are 7% for all beam energies, dominated by the uncertainty 
on the ratio of proton elastic to proton dissociative cross sections. This is added in quadrature 
to the uncertainty of 3(4)% on the luminosity measurement to obtain the total normalisation 
uncertainty of 7.6(8.1)% for the 920 GeV (LME) data. 

A full decomposition of the systematic errors on the measured cross sections is given in 
tables [21 [3] and m Correlated sources of uncertainty that are always smaller than 2% and are 
never the dominant correlated source in a single bin are omitted. For the LME data, the precision 
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of the cross section measurements is statistically limited in the region of greatest sensitivity to 
F£ at high y. Elsewhere in the LME data, the systematic errors are of similar size to the 
statistical errors. The precision of 4% reached in the best-measured regions for the 920 GeV 
data is the highest accuracy achieved in HI measurements of a^ to date. The 920 GeV data are 
limited by the systematic uncertainties throughout the measured range, the dominant source of 
systematic uncertainty varying with the kinematics. The largest correlated uncertainty at low 
xjp comes from the modelling of the LAr noise, with the vector meson simulation also playing 
an important role. At low (3 (high y), where Fj^ is measured, the largest sources of uncertainty 
are the photoproduction background subtraction, the efficiency of the riog cut and the model 
dependence arising from the F£ treatment in the MC simulation. The uncertainty arising from 
imperfect knowledge of the bin-centre corrections can also be large, typically at large /3, low 
Xjp or low Q'^. 



4.7 Extraction of Ff 

The separation of F2 and F£ follows a similar procedure to that which was used to extract their 
inclusive counterparts F2 and Fl [TS]. The diffractive reduced cross section is integrated over 
the My and t ranges given in equation [51 The uncertainty on correcting an individual dataset 
to that range is large (7%) but strongly correlated between datasets. The residual difference in 
normalisation between the three datasets after all corrections is determined from comparisons of 
a^ at low y to be 2%. In order to extract F^ optimally, the cross sections are normalised to the 
HI 2006 DPDF Fit B result in a range where the sensitivity to F^ is minimal, but the statistical 
precision and kinematic overlap of the data is still sufficient. Data in the range Q^ > 7 GeV^, 
xp = 0.003 and y < 0.38 (0.3 and 0.3) for the 460 (575 and 920 GeV) datasets are used, 
yielding normalisation factors of 0.97, 0.99 and 0.97, respectively. As the published data at 
820 GeV were included in the analysis of the data used as input to the HI 2006 DPDF Fit B, 
they are already consistently normalised. 

Following this normalisation procedure, the diffractive longitudinal structure function Fj^ 
can be extracted directly from the slope of cr^ as a function of y'^/Y^ for each set of Q^, xp 
and /3 values. A linear fit is performed, taking only the statistical errors, dgtat, into account in 
order to calculate the statistical uncertainty on Fj^ . The fit is repeated, adding the statistical and 
uncorrected errors in quadrature, 5stat+unc, to calculate the measured value of F^ and the sum 
of its combined statistical and uncorrected errors. For each correlated systematic error source, 
each of the cross section points is adjusted according to the positive and negative shifty and 
the fit is repeated using 5stat+unc for the errors on the cross section points. The error on Fj^ is 
taken as half of the difference between fits to the positive and negative shifted data points. All 
of these correlated errors are added in quadrature with 5stat+unc to give the total error on F^ . 
The normalisation uncertainty on the value of F^ is set by the normalisation uncertainty on the 
cross section measurements and is therefore 8.1%. 

As only bin-averaged cross sections are available at the highest /3 > 0.9, Ff is not extracted 
in that region. 



^ In fits which include the published 820 GeV data, a more conservative approach is used whereby the 820 GeV 
data remain fixed. This results in a larger variation in the slope of a^ as a function of y^ j Yj^ with a correspondingly 
larger uncertainty on F^ . 
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4.8 Extraction of R^ and the ratio R^/R 

The photoabsorption ratio for diffraction, R^ = Fj^ /{F^ — F^), is extracted from linear fits to 
the data by reparametrising equation|4]such that R^ and (F^ — F^) become the free parameters 
of the fit: 

<y? = {F^ - fE) + R"" ■ {fE - fE) ■ (1 - yVy+)- (11) 

The error on R^ is calculated in the same way as for FE , detailed in section 14771 The normali- 
sation uncertainty cancels in this ratio. 

In order to calculate the ratio of i?^ to its inclusive counterpart R = Fl/{F2~Fl), the value 
of R is extracted from the present data using a similar procedure to that used for R^ described 
above. Only data with Q'^ > 7 GeV^ are used, where inclusive measurements are made at all 
beam energies in this analysis. The statistical correlations between the inclusive and diffractive 
measurements are neglected and the systematic errors are assumed to be dominated by the error 
on R^. Similarly to R^, there is no normalisation uncertainty on the ratio R^ / R. 



5 Results 

The measured diffractive reduced cross section values and their errors are given in tables [2l [3] 
andlH Figure [5] shows the reduced cross section as a function of (3 at fixed xp and Q^ for the 
LME, 820 GeV and 920 GeV datasets. Also shown is the prediction of HI 2006 DPDF Fit B, 
which in general describes the data well at Q'^ > 11.5 GeV^. Deviations of the measured cross 
sections from the F^ predictions at low /3 are evident in the LME data, where the highest y 
values are accessed, notably at Q^ = 11.5 GeV^ and xp = 0.003. This shows the sensitivity 
of the LME data to Ff . The extrapolation to lower Q^ of HI 2006 DPDF Fit B, which only 
included data with Q^ > 8.5 GeV^, is also compared with the Q"^ = A GeV^ data. The fit is 
known to significantly undershoot the published 820 GeV data in this region [|3]|, an observation 
which is reproduced for the new measurements. 

The new data at xjp = 0.0005, Q^ = 11.5 GeV^ and xjp = 0.003, Q^ = 44 GeV^ include 
the highest /3 measurements obtained by HI to date. They are in remarkably good agreement 
with the extrapolation of HI 2006 DPDF Fit B and support the hypothesis that a^ — ;■ as 
/3 — )■ 1. There is thus no evidence in this region for a large higher twist FE contribution 
[1211261123. 

The extraction of Ff via linear fits to the y'^/Y^ dependence of the reduced cross section at 
different beam energies and fixed Q^, (3 and xp is shown in figure l6l The largest lever arm in 
y'^ /Y+, and therefore the highest sensitivity to FE , is at the lowest (3. The data are consistent 
with a linear dependence of a^ on y^ jY^, with a significant tendency for o^ to decrease as 
y'^/Y^ increases for most Q^, xp and /3 values. The values of FE and their errors are given in 
table 121 

The measurements of Ff , at fixed values of Q^ and xp, are shown as a function of (3 in fig- 
urelTl Significantly non-zero measurements of FE are made for all values of Q^ and xp and five 
FE points are greater than zero by more than 3a. The data are compared with the predictions 
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of the HI DPDF Fits A and B JSl and with the Golec-Biernat & Luszczak model |l29l| (section 
[3]). All three models are consistent with the data, although there is a tendency for the measure- 
ments to lie above the predictions. Although the prediction of [29| lies significantly above both 
Fit A and Fit B at large /3, the experimental precision is insufficient to distinguish between the 
models. The measured values of F^ are also shown in figure |71 The F^ measurements agree 
well with the predictions of HI DPDF Fit B for Q^ > 11.5 GeV^. Within the uncertainties, all 
measurements are consistent with the hypothesis that < F^ < F^. 

A summary of the F£ measurements is given in figure [8l where the data points from all 
five Q^ and xp values are shown as a function of (3 and compared with the HI 2006 DPDF 
Fit B prediction. In order to remove the significant dependence on xp, the F^ points have been 
divided by a factor fip/p, taken from [3j, which expresses the measured xp dependence of the 
data, assuming proton vertex factorisation. The remaining discontinuities in the prediction are 
due to its Q^ dependence. The F^ data cover a large range in longitudinal fractional momentum 
0.033 < (3 < 0.7 and are compatible with the predicted slow decrease with increasing /3. 
The data have a tendency to lie above the prediction although the precision is limited. The 
most significantly positive F^ measurements lie in the region (3 < 0.5, which contrasts with 
models of diffraction such as [24,26], which do not include leading twist contributions from 
longitudinally polarised photons. 

The measurement of R^ is shown as a function of (3 in figure [9l Data with \R^\ > 50 and a 
relative uncertainty larger than 100% are not shown. The data are compatible with the prediction 
based on HI 2006 DPDF Fit B, though they are also consistent with other models. The data 
at Q^ = 11.5 GeV^ indicate that the longitudinally and transversely polarised photon cross 
sections are of the same order of magnitude (R^ ~ 1 and F^ ~ 2F^). At Q'^ = 44 GeV^, 
where larger (3 values are accessed, there is a tendency for the data to lie above the prediction, 
which tends to zero as (3 — t- 1 . There is no evidence for the steep rise in R^ which might be 
expected at large (3 if configurations similar to vector meson electroproduction were dominant 
in this region. The values of R^ and their errors are given in table [6l 

The relative importance in inclusive and diffractive scattering of the longitudinally po- 
larised photon cross section compared with its transverse counterpart is investigated via the 
ratio R^ /R, shown as a function of x in figure \Wi Only data with Q"^ > 7 GeV^, where a 
measurement of R is possible in this analysis, are used. Data with \R^ /R\ > 20 and a relative 
uncertainty greater than 100% are not shown. The ratio data suggest that the longitudinally 
polarised photon contribution plays a larger role in the diffractive than the inclusive case. Aver- 
aged over all data, R^ /R = 2.8 ± 1.1. The data are well reproduced by the ratio of predictions 
from HI 2006 DPDF Fit B and an HI fit to inclusive DIS data, HI PDF 2009 [46]. At high Q^ 
corresponding to high x and therefore (3, the prediction decreases towards zero as x — )■ 1. The 
data are consistent with such a decrease with increasing /3 within large experimental uncertain- 
ties. 



6 Conclusions 

First measurements of the diffractive reduced cross section at centre-of-mass energies y/s of 
225 and 252 GeV are presented, together with a precise measurement at ^/s of 319 GeV. The 
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reduced cross section is measured in the range of photon virtualities 4.0 < Q^ < 44.0 GeV^ 
and of the longitudinal momentum fraction of the diffractive exchange 5 • 10^"^ < xjp < 3 • 10^^. 
The reduced cross section measurements agree well with predictions derived from leading twist 
NLO QCD fits to previous HI data throughout the kinematic range. The data at high and 
medium inelasticity y are used to extract the first measurement of the longitudinal diffractive 
structure function F^ . There is a tendency for the predictions to lie below the F^ data, but the 
data are compatible with HI 2006 DPDF Fit A and Fit B as well as with a model which includes 
a higher twist contribution at high (3, based on a colour dipole approach. The procedure also 
allows a simultaneous extraction of F2 , independently of assumptions made on F^, in the 
same kinematic range. The F2 measurements agree well with the predictions of HI DPDF 
Fit B for Q^ > 11.5 GeV^. Within the uncertainties, all measurements are consistent with the 
expectation that Q < F^ < F^ . 

The ratio i?^ of diffractive cross sections for longitudinally to transversely polarised pho- 
tons is measured in the same kinematic range as F^ . At fixed Q^ and xp, this ratio is relatively 
flat as a function of /3 and suggests that the cross sections for the two polarisation states of the 
photon are of comparable size. The ratio of R^ to its inclusive scattering counterpart, R, is 
extracted in the region Q^ > 11.5 GeV^. The R^ / R data indicate that the longitudinally po- 
larised photon cross section plays a larger role in the diffractive than in the inclusive case. The 
R^ and R^ /R measurements are well reproduced by the predictions based on HI 2006 DPDF 
Fit B and the HI PDF 2009 inclusive PDF set. 
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Table 2: The diffractive reduced cross section af at -y/s = 225 GeV, multiplied by xp, measured with the 460 GeV data, at fixed values 
of xjp, Q^ and (5. At the largest /3, the bin-averaged cross section is given together with the lower and upper bin boundaries. The statistical 
i^stat), uncorrected (5„„c) and sum of all correlated {5cor) uncertainties are given together with the total uncertainty {5tot)- The other columns 
show the individual correlated uncertainties, which are due to the positron energy scale {5eie), the positron polar angle measurement {be), the 
LAr noise subtraction {5noi), the hadronic SpaCal energy scale {5 spa), the efficiency of the logarithmic energy- weighted cluster radius cut 
((5r,^ ), the charge asymmetry of the photoproduction background {5asy), the model uncertainty due to the influence of Ff {5mod), the model 
uncertainty on the underlying (5 and xjp distributions {5p, d^jp), the influence of resonant (5^m) and QED Compton {5com) contributions and 
finally the parametrisation choice for the bin centre corrections {Sbcc)- A minus sign indicates that a source is anti-correlated with a change 
in the cross section. All uncertainties are are given in per cent. The normalisation uncertainty of 8.1% is not included. 
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Table 3: The diffractive reduced cross section a^ at ^/s = 252 GeV, multiplied by xjp, measured with the 575 GeV data, at fixed values 
of X]p, Q"^ and (3. At the largest (3, the bin-averaged cross section is given together with the lower and upper bin boundaries. The statistical 
{5stat), uncorrelated (Sunc) and sum of all correlated (5cor) uncertainties are given together with the total uncertainty i6tot)- The other columns 
show the individual correlated uncertainties, which are due to the positron energy scale {6eie), the positron polar angle measurement {6g), the 
LAr noise subtraction {6noi), the hadronic SpaCal energy scale {6spa), the efficiency of the logarithmic energy-weighted cluster radius cut 
(5r,^ ), the charge asymmetry of the photoproduction background (Sasy), the model uncertainty due to the influence of Ff (Smod), the model 
uncertainty on the underlying (3 and xp distributions (5^, 5xjp), the influence of resonant (5vm) and QED Compton (Scom) contributions and 
finally the parametrisation choice for the bin centre corrections (Sbcc)- A minus sign indicates that a source is anti-correlated with a change 
in the cross section. All uncertainties are given in per cent. The normalisation uncertainty of 8.1% is not included. 
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Table 4: The diffractive reduced cross section a^ at y/s = 319 GeV, multiplied by xp, measured with the 920 GeV data, at fixed values 
of X]p, Q'^ and (3. At the largest (3, the bin-averaged cross section is given together with the lower and upper bin boundaries. The statistical 
{5stat), uncorrected {5unc) and sum of all correlated {5car) uncertainties are given together with the total uncertainty {6tot)- The other columns 
show the individual correlated uncertainties, which are due to the positron energy scale (^e^e), the positron polar angle measurement {6g), the 
LAr noise subtraction (5„oj), the hadronic SpaCal energy scale i5spa), the efficiency of the logarithmic energy-weighted cluster radius cut 
{6ri^ ), the charge asymmetry of the photoproduction background {6asy), the model uncertainty due to the influence of Ff {5rnod), the model 
uncertainty on the underlying (3 and xjp distributions (5^, Sxp), the influence of resonant (6^^) and QED Compton (Scom) contributions and 
finally the parametrisation choice for the bin centre corrections {5hcc)- A minus sign indicates that a source is anti-correlated with a change 
in the cross section. All uncertainties are given in per cent. The normalisation uncertainty of 7.6% is not included. 
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Table 5: The diffractive structure functions Ff and F^ multiplied by X]p, at fixed values of xjp, Q^ and (3. The statistical uncertainty {5stat), 
the sum of the statistical and uncorrelated uncertainties {5stat+unc) and the sum of all correlated uncertainties {6cor) are given together with 
the total uncertainty {5tot)- Absolute uncertainties are given. The normalisation uncertainty of 8.1% is not included. 
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Table 6: The ratio Rd of the cross sections for longitudinally to transversely polarised photon cross sections, at fixed values of xp, Q^ and 
l3. The sum of the statistical and uncorrected uncertainties (6stat+unc) and the sum of all correlated uncertainties (5cor) are given together 
with the total uncertainty {6tot)- Absolute uncertainties are give. 
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Figure 2: The energy distributions of the scattered positron candidates for the 460 GeV (left) 
and 575 GeV (right) data. The data shown as points are compared with the sum of the diffrective 
DIS MC simulation and background estimates (open histogram). The light-filled histogram 
shows the photoproduction background estimate from data, the dark-filled histogram is the sum 
of the QED Compton and inclusive DIS backgrounds, taken from MC simulations. 
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Figure 3: The quantity Sj(-E — pz)i summed over all final state particles for the 460 GeV (left) 
and 575 GeV (right) data at high y. The data after background subtraction are shown as points, 
compared with the MC simulation shown as a histogram. The shaded area shows the effect of a 
variation of the hadronic SpaCal energy scale by its uncertainty of 5%. 
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Figure 4: Distributions of the kinematic quantities y (top), (3 (middle) and log(a;jp) (bottom) 
for the 460 GeV (left), 575 GeV (middle) and 920 GeV (right) datasets. The data are shown 
as points compared with the sum of the MC simulation and background estimates (open his- 
togram). The light-filled histogram shows the photoproduction background estimate from data, 
the dark-filled histogram is the sum of the QED Compton and inclusive DIS backgrounds, ob- 
tained from MC simulations. 
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Figure 5: The diffractive reduced cross section a^ multiplied by xp as a function of (3 at fixed 

g2 and xjp for (from left to right) the 460 GeV, 575 GeV, 820 GeV and 920 GeV datasets. The 

data are compared with the predictions of HI 2006 DPDF Fit B (solid line), which is indicated 

as dotted beyond the range of validity of the fit. The dashed and dashed-dotted lines represent 

the contribution of F^ , which is the same for each beam energy. The inner error bars represent 

the statistical errors on the measurement, the outer error bars represent the statistical and total 

systematic uncertainties added in quadrature. The normalisation uncertainties of 7.6(8.1)% for 

the 920(460, 575) GeV data are not shown. 
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Figure 6: The diffractive reduced cross section a^ multiplied by xp as a function of y'^/Y^ at 
fixed Q^, xp and (5. The inner error bars represent the statistical uncertainties on the measure- 
ment, the outer error bars represent the statistical and total systematic uncertainties added in 
quadrature. The normalisation uncertainty is not shown. Up to four beam energies are shown, 
where the lowest |/^/l+ point is given by the 820 GeV data for Q^ = A GeV^ and by the 
920 GeV data at higher Q"^. The linear fits to the data are also shown as a solid line, the slope of 
which gives the value of F^ . The predictions and extrapolated predictions of HI 2006 DPDF 
Fit B are shown as dashed and dotted lines, respectively. 
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Figure 7: The diffract! ve structure functions F^ and F^ multiplied by xjp as a function of (3 
at fixed Q"^ and xjp. The F^ data are shown as filled points, compared with the predictions 
of HI 2006 DPDF Fit A (dashed line), Fit B (solid line) and the Golec-Biernat and Luszczak 
model (dashed and dotted line). The measurements of F,p (open points) are compared with 
the prediction of HI 2006 DPDF Fit B (long dashed line). The inner error bars represent the 
statistical uncertainties on the measurement, the outer error bars represent the statistical and 
total systematic uncertainties added in quadrature. The normalisation uncertainty of 8. 1% is not 
shown. 
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Figure 8: The diffract! ve longitudinal structure function Ff, divided by a parametrisation of 
the xjp dependence of the reduced cross section fjp/p [[3l, as a function of /3 at the indicated 
values of Q"^ and xp. The data are compared with the predictions of HI 2006 DPDF Fit B 
(red line), which is indicated as dashed beyond the range of validity of the fit. The inner error 
bars represent the statistical uncertainties on the measurement, the outer error bars represent the 
statistical and total systematic uncertainties added in quadrature. The normalisation uncertainty 
of 8.1% is not shown. 
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Figure 9: The ratio R^ of cross sections for longitudinally to transversely polarised photons, 
as a function of f3 at the indicated values of xp and Q^. The data are compared with the 
predictions of HI 2006 DPDF Fit B, indicated as dashed beyond the range of validity of the fit. 
The inner error bars represent the statistical uncertainties on the measurement, the outer error 
bars represent the statistical and total systematic uncertainties added in quadrature. 
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Figure 10: The ratio of R^ /R as a function of x at the indicated values of Q'^ and xp. The data 
are compared with the predicted ratio using HI 2006 DPDF Fit B / HI PDF 2009 (solid line). 
The error bars represent the statistical and systematic uncertainties added in quadrature. 
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